Head shock vs Mach cone: azimuthal correlations from 2— )-3 parton processes in 

relativistic heavy-ion collisions 
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We study the energy-momentum deposited by fast moving partons within a medium using hn- 
earized viscous hydrodynamics. The particle distribution produced by this energy-momentum is 
computed using the Cooper-Frye formalism. We show that for the conditions arising in heavy-ion 
collisions, energy momentum is preferentially deposited along the head shock of the fast moving 
partons. We also show that the double hump in the away-side of azimuthal correlations can be pro- 
duced by two (instead of one) away-side partons that deposit their energy-momentum along their 
directions of motion. These partons are originated in the in-medium hard scattering in 2 — > 3 pro- 
cesses. We compare the results of the analysis to azimuthal angular correlations from PHENIX and 
show that the calculation reproduces the data systematics of a decreasing away-side correlation when 
the momentum of the associated hadron becomes closer to the momentum of the leading hadron. 
This scenario seems to avoid the shortcomings of the Mach cone as the origin of the double-hump 
structure in the away-side. 

PACS numbers: 25.75.-q, 25.75.Gz, 24.10.Lx, 24.10. Nz 



I. INTRODUCTION 

Azimuthal angular correlations have provided a povifer- 
ful testing ground to elucidate the propagation properties 
of fast partons within the medium created in high-energy 
heavy- ion reactions 0,0 ■ The main features of these cor- 
relations can be summarized as follows: when the lead- 
ing and the the away-side particles have similar momenta, 
the correlation shows a suppression of the away-side peak, 
compared to proton collisions at the same energies. How- 
ever, when the momentum difference between leading 
and away-side particles increases, either a double peak or 
a broadening of the away-side peak appears, while nei- 
ther of these are present in proton collisions at the same 
energies [|. The peak region is known as the shoulder 
and the region between the peaks is known as the head. 
For low momentum particles, the jet multiplicity is larger 
in the shoulder than in the head region [2] . 

Since the structures in the away-side are best seen 
for low momentum particles, there were explanations 
based on the emission of sound modes caused by one 
fast moving parton 0, Q, the so called Mach cones. 
Nevertheless, it has been argued that such interpreta- 
tion is fragile, since the jet-medium interaction produces 
also a wake whose contribution cannot be ignored 0] . 
Furthermore, it has recently been shown that it is un- 
likely that the propagation of one high energetic particle 
through the medium leads to a double peak structure in 
the azimuthal correlation in a system of the size and fi- 
nite viscosity relevant for heavy ion collisions, since the 
energy-momentum deposition in the head shock region is 
strongly forward-peaked Q]. Moreover, by using a real- 
istic multiple gluon emission for the parton shower pro- 
duced by an in-medium moving parton, the overlapping 
perturbations in very different spatial directions wipe out 



any distinct Mach cone structure Q. 

The more widely accepted explanation for the double 
peak/broadening of the away-side is given nowadays in 
terms of initial state fluctuations of the matter density in 
the colliding nuclei. These fluctuations are then shown to 
give rise to an anisotropic flow of partially equilibrated, 
low momentum particles, within the bulk medium. Hy- 
drodynamic descriptions of this scenario have for instance 
successfully reproduced the experimental (lo| . How- 
ever, it has also been shown that there is a strong connec- 
tion between the observed away-side structures and the 
medium's path length llj. The connection is expressed 
through the dependence of the azimuthal correlation on 
the trigger particle direction with respect to the event 
plane in such a way that, for selected trigger and associ- 
ated particle momenta, the double peak is present/absent 
for out of plane/in plane trigger particle direction. A fi- 
nal state effect rather than an an initial state one seems 
more consistent with this observation. 

The medium's response to a fast moving parton can 
be described by the energy and momentum that the par- 
ton deposits. This energy-momentum is then converted 
into particles upon hadronization. When the interaction 
produces low momentum particles, a hydrodynamical de- 
scription of the medium's response seems appropriate. 
In this framework, it is natural to explore the types of 
medium excitations generated by fast moving partons 
and the way these in turn are responsible for the shape 
of the away-side in the angular correlations. When the 
sound (that is the longitudinal modes ) are prominently 
excited, the double peak in the away-side could emerge 
as a Mach cone structure, that is, as energy-momentum 
mainly deposited on the sides of the path travelled by 
the parton. This is illustrated in the left panel of fig. [T] 
However, when the modes that get prominently excited 




FIG. 1: (Color on-line) Schematic representation of the possible origin of the away-side double-hump. The left panel illustrates 
the common picture whereby one away-side parton deposits energy-momentum, exciting the sound modes and producing a 
Mach cone. The right panel illustrates our approach whereby two away-side partons generate each a head-shock and deposit 
energy-momentum preferentially along their directions of motion. 



are the wake, (that is the transverse modes), energy and 
momentum is preferentially deposited along the direction 
of the travelling parton. If a double peak is to emerge, 
the alternative scenario is that the structure is produced 
by two travelling partons in the away-side. This is illus- 
trated in the right panel of fig. [TJ 

Two partons in the away-side are produced in processes 
where two initial state partons scatter into three final 
state ones, the so called 2 — >■ 3 processes. In this case the 
cross section is smaller than the one for 2 —s- 2 processes, 
at least by one power of a^. However, consider a scat- 
tering event that deposits a given amount of energy in 
the away-side. If such event comes from a 2 — >■ 2 process, 
the away-side parton will have a larger energy than each 
of the two partons, when these last come from a 2 — > 3 
process. Since the parton distribution is a fast falling 
function of the parton energy, the extra power of is 
partially compensated by the larger abundance of par- 
tons with smaller energy. This is illustrated in fig. [21 For 
2—^3 processes, when the dominant energy-momentum 
deposition happens through transverse modes, conserva- 
tion of momentum at the parton level gives rise to a dis- 
tinctive angular dependence in the azimuthal correlation 
whereby, the angular difference between the peaks in the 
away-side is close to 27r/3 radians. 

The computation of hadron events produced by 2 — >■ 3 
parton processes in the context of azimuthal correla- 
tion functions has been put forward and explored in 
Refs. [13, US], using the leading order QCD matrix el- 
ements. Those studies were made for fragmentation out- 
side the medium. In this work we compute the hadron 
multiplicity produced by two partons in the away-side 
coming from 2^-3 processes for the case where all its 
energy-momentum is deposited into the medium. We 
use linearized viscous hydrodynamics to compute the 



medium's response. The work is organized as follows: In 
Sec.|II]we first review the general framework to compute 
the particle distribution stemming from a given energy- 
momentum deposited by a fast moving parton within the 
medium, using linearized viscous hydrodynamics. This 
energy-momentum is then converted into a parton multi- 
plicity using the Cooper- Frye formalism . We discuss 
the different shapes for the azimuthal correlations that 




FIG. 2: Schematic representation of how a single away-side 
parton is less abundantly produced than two away-side par- 
tons carrying the same total momentum, given that the mo- 
mentum distribution function rapidly decays with momen- 
tum. 
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are obtained when varying the strength of the wake and 
sound mode contributions. In Sec. IIIII we solve the hy- 
drodynamic equations for the transverse and longitudinal 
modes. We work in the limit where the patron's veloc- 
ity is larger than the sound velocity. In Sec. |TV] we use 
these solutions to compute the azimuthal correlations by 
convoluting the pQCD probability to produce two away- 
side partons in 2 — > 3 processes with the multiplicity 
obtained from the excess of energy-momentum produced 
by the moving partons followed by fragmentation. We 
finally summarize and conclude in Sec. |Vl 



II. PARTICLE DISTRIBUTION 

We compute the particle's multiplicity as given by the 
Cooper-Frye formula 



E 



dN 



1 



dj:^p^[f{p-u)-f{po)], 



(1) 



where f{p ■ u) — f{po) is the phase-space disturbance 
produced by the fast moving parton on top of the equi- 
librium distribution /(po): with and represent- 
ing the freeze-out hyper-surface and the particle's mo- 
mentum, respectively. The medium's total four-velocity 
m'' = Uq + (5u^ is made out of two parts: the background 
four- velocity Uq and the disturbance (5u''. This last con- 
tribution is produced by the fast moving parton and can 
be computed using viscous linear hydrodynamics once 
the source, representing the parton, is specified. For a 
static background (which hereby wc assume) and in the 
linear approximation, can be written as 



1, 



eo(l + c2) 



(2) 



where the spatial part of the medium's four-velocity, 
u = g/eo(l + Cg), is written for convenience in terms of 
the momentum density g associated to the disturbance, 
with eo and Cs the static background's energy density and 
sound velocity, respectively. We focus on events at cen- 
tral rapidity, y ~ and take the direction of motion of 
the fast parton to be the z axis and the beam axis to be 
the X axis. With this geometry, the transverse plane is 
the y — z plane and therefore, the momentum four- vector 
for a (massless) particle is explicitly given by 



= {E-,Px,Py,Pz) 

= (pt, 0, Pt sin 4>,pt cos ( 



(3) 



where is the angle that the momentum vector p makes 
with the z axis. We use Bjorken's geometry, thus 

d'^p = pxdpTd(f>dpx 

Px = PTsinhy 

dpx = Pt coshy dy 

E = pxcoshy, (4) 



and therefore 
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dN 



dN 



d^p pxdpT d(j)dy 



(5) 



For simplicity we consider a freeze-out hyper-surface of 
constant time 



d^. 



(dV, 0,0,0). 



(6) 



Therefore, using Eqs. ([T]) and ([S|), the particle azimuthal 
distribution around the direction of motion of a fast mov- 
ing parton within the px interval p^'" < pr pf""" 
given by 



dN 
dy 



y~0 



ArAy 
"(2^ 



Pt 



dpTPr 



X d\[f{p-u)-f{po)], 



(7) 



with At the freeze-out time interval and where we have 
assumed perfect correlation between the space-time ra- 
pidity 77 and y to substitute A77 by Ay. We assume that 
the equilibrium distribution is of the Boltzmann type. In 
this way, we have 



fiPo) 
f{p ■ u) 



exp[pT/To] 
Pt 



exp 
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gy sm ( 



gz cos ( 



(8) 



where Tq is the background medium's temperature and 
6T is the change in temperature caused by the passing of 
the fast parton. Assuming that the energy density and 
temperature are related through Boltzmann's law 



eocT\ 



one gets 



ST 



5e_ 
4eo' 



(9) 



(10) 



Since for the validity of linearized hydrodynamics, both 
5e and g need to be small quantities compared to Eq, we 
can expand the difference f{p ■ u) — f{pa) to linear order. 
Using Eqs. (g]) and dTO]) we get 



f{p ■ u) - fipo) 



Pt_ 
To 

il 
4eo 



exp [~pt/Tq] 

gy sm(j) + gz cos( 
eo(l + c2) 



(11) 



therefore the particle azimuthal distribution around the 
direction of motion of a fast moving parton is given, in 
the linear approximation, by 



dN 
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Notice that the shape of the distribution around the di- 
rection of the fast moving parton depends on the quan- 
tities 



Gy 



(13) 



When Gy > Gz the distribution is dominated by the 
sint/) factor, giving rise to two peaks away from (f> = 0. 
However, when Gz > Gy the opposite happens and the 
distribution is dominated by the cos factor with the two 
peaks close to cj) ~ 0. These peaks become a single peak 
in the extreme case where Gz ^ Gy. This is illustrated 
in fig. [3] 

We now proceed to show that when the velocity of the 
moving parton is larger than the speed of sound, Gy is 
mostly made out of longitudinal (sound) modes, whereas 
Gz is mostly made out of wake (transverse) modes and 
that the latter dominates the former, giving rise to a par- 
ticle distribution centered around the direction of motion 
of the moving parton. 



III. LINEAR VISCOUS HYDRODYNAMICS 

In order to compute the gy and gz components of the 
momentum density vector g which is deposited into the 
medium by the moving fast parton, we resort to using lin- 
earized viscous hydrodynamics. Assuming that the dis- 
turbance introduced by the parton is small, the medium's 
energy-momentum tensor can be written as 



(14) 



where 60'^'^ is the perturbation generated by the parton 
and Qq'^ is the equilibrium energy-momentum tensor of 
the underlying medium. Therefore each of these compo- 
nents satisfy the equations 



0, 



(15) 



where J'^ represents the source of the disturbance, which 
in this case corresponds to the fast moving parton. Here 
we consider that the parton can be represented by a lo- 
calized disturbance of the form 



(16) 



where {dE/dx) is the energy loss per unit length and 

i>''=7(l,v), (17) 

with 7 the Lorentz gamma factor. Effects of a finite 
source structure have been studied in Ref. [l5| where it 
is found that differences in the energy density deposi- 
tion between a finite extent source and a localized one 
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FIG. 3: Schematic representation of the particle distribution, 
arising from the Cooper-Frye formula. The distribution is 
centered around the direction of motion {(j) — 0) of the parton 
that deposits the energy-momentum in the medium. Since 
the distribution is a sum of a sm{(j)) and a cos((j!)), its shape 
depends on the relative strength of their respective coefficients 
Gy and Gz- From top to bottom, we show first the sm{4)) and 
a cos{(f>) functions, the case where Gy > Gz, the case where 
Gz > Gy and finally the case where Gz ^ Gy. 
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exist only close to the source. Since a hydrodynamical 
description is anyway valid only for large distances from 
the source, compared to the transport mean free path, 
here we will consider that the description of a source as 
localized exactly at the position of the fast parton suf- 
fices. 

To describe the propagation of the disturbance caused 
by the source, Eqs. (jisp are solved by considering that 
the energy-momentum tensor 8^"^ consists of a piece that 
describes an isotropic fluid 



(18) 



and a disturbance 60'^'^ caused by the source, whose ex- 
plicit components to first order in shear (rj) and bulk ((^) 
viscosity, are given by 

60°" = 6e, 

- Cd^.V.g, (19) 

where e(t, x) = eo + 6e{t, x) with eg corresponding to the 
energy density of the background fluid, Se to the energy 
density associated to the disturbance and 



4?/ 



3eo(l + c2 



(20) 



is the sound attenuation length. 

In the linear approximation and vanishing bulk viscos- 
ity (l5| , the dynamical description of the propagation of 
the disturbance is given by the first of Eqs. (|15p . whose 
explicit components can be written as 



(21) 



by [If 



(5e(k,a;) 
gL(k,w) 
gT(k,cj) 



ik ■ J(k,w) + J0(k,w)(iw - r,fc2) 
z[|^k.J(k,c^)+c^J°(k,c^)]k 



= g - gL = 



iJT(k, uj) 



,(24) 
(25) 
(26) 



where we have written the momentum density vector in 
terms of its transverse and longitudinal components with 
respect to the Fourier mode k, namely 



gL + gT 



(27) 



with the definition of longitudinal and transverse compo- 
nents of any vector cr given by 



(^k 

cr - (Tl. 



(28) 



Note that the source term J'^ in Eq. ([T5)) is Fourier trans- 
formed to 



J''(k,a;) = (27r)7 



/dE 
\ dx 



(5(k • V — uj)v'^ , 



(29) 



where we have assumed a constant energy loss per unit 
length. 

From Eqs. ([Mjl - ipS]) one can obtain the space time 
solutions for 5e(x, t) and g(x, t) upon use of Eq. (|22p . 
The corresponding expressions are 



. ^ ,^ ^ ,'dE\ r d^k 

gT(x,t) = Z7(2^) I — I 



dw 



dx J J (27r)3 J (27r) 



AXL-yi—iujt 



(k - J)k 



(5(k • V — w) 

UJ 4- ijTsk'^ 



(30) 



These equations can be readily solved in momentum 
space. We define the Fourier transform pair /(x, t) and 
/(k,w) as 

= (2^ / '^'^ / (22) 



Using Eq. (|22]) into Eqs. (|2T]), together with Eqs. (fT9| . 
we obtain 

— za;(5e + ik • g = J*^, 
-^^g» + ^c2fc»5e+ Jr,(fcV + y(k-g)) = J^ (23) 

Thus the perturbed energy density, 5e and momentum 
density components g*, can be obtained by solving the 
algebraic system of Eqs. (|23p and the solutions are given 



/ X N fdE\ f d^k 

- *^(2^)Ujy (2^ 



duj 



'i\i.--x.—iujt 



X k 



(5(k • V — w) 



(jj- — c%k^ + iTsLok" 



(31) 



and 



5e{yL,t) = i7(27r) 



dE 



d^k 



duj 



(27r)3 J (27r) ' 



Lj^ — c?gk'^ + iT sUik"^ 

(32) 



[ik • V + zcli — r^fc^] 



In order to compute the integrals in Eqs. (p0 l) - ([32|) we 
use cylindrical coordinates with the fc^ direction along 
the direction of motion v of the fast parton. It is easier 
to start with the components of gT- After carrying out 
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the frequency and angular integration, we get for the z 
component 



dk 



{z — Vt) 



dxj Jo (2^)2 (27r) 

kpojkTXT) , . 

(fcf + fc2)[t;A:,+zfr«(fc2 + fc2)]' ^ ) 

where Jq is a Bessel function and xt = \p\P' is the 
distance from the parton along the transverse direction 
(along the y axis in the geometry we are using). The in- 
tegration over k^ is performed using contour integration. 
For causal motion {z — vt > 0) we close the contour on 
the lower half fcz -plane. The poles that contribute are 
located at 



In a similar fashion we get 



V47r 



dE 
dx 



2v 

3f: 



(^2 +^2)2 

d.s^s{s + 2)Ji (/3 + 2)) e-"' 



2v 



The computation of the components of is more in- 
volved due to the analytic structure of the poles in the 
complex /cz-plane. Let us fist compute the z component. 
After carrying out the frequency and angular integration 
we get 



-iki 



3r, 



(34) 



After carrying out the k^ integration, the remaining inte- 
gral can be expressed in terms of the dimensionlcss quan- 
tities 



1 i. 



1^-^*1/ (If 



/ /^3r. 

XT — 

\ 2v 



(35) 



as 



(St)^ = 



= -CL 



(36) 



The first term in Eq. ([36]) can be analytically integrated. 
For numerical purposes, it is more convenient to rewrite 
the second term after the change of variable 



(37) 



The final result is 
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(a2 + p^) 
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(2^)2 (27r) 



£_ ifcs (z — vt) 



kTJo{kTXT)k:,[k'^ + ^{kl + fc|)] 

(fc2 + fc2)[fc2 + (,IJ.fc^_£|)(fc2+fc2)]- 



(40) 



Notice that the integrand contains the parameter 
Cg/v'^. For a fast moving (masslcss) parton u ~ 1 and for 
a relativistic gas, Cs — ■\/l/3. Therefore, for conditions 
close to the ones after a heavy- ion reaction c^/v^ < 1 and 
we can expand the integrand in Eq. (|40p in this parame- 
ter. To first order in c^/w^, we get 

, , ■ ,n \ ( dE\ f°° dkx , s 
(g.). = ^7(2^) (^^j I j^krMkrxr) 



dkz 

W) 



y(j,2gifc- {z—vt) 



[kl + kl)% + {^^){kl + kl)] 



c: ^ 2kS'-t){ki + k'T)^ 

2 



V 



ik- (z — vt) 



[kz + i^^m + kl)] 



(41) 



For causal motion {z — vt > 0) we close the contour on 
the lower half fc^-plane. The poles that contribute are 
located at 



-ikx 



(42) 



The remaining integral over fc^ can be expressed in 
terms of the dimensionlcss quantities defined in Eq. (|35p . 
For numerical purposes, it is more convenient to rewrite 
this integral after the change of variable in Eq. ([57]) . The 
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FIG. 4: (Color on-line) Three dimensional plots and the corresponding contour plots for (top to bottom) /g^ _ , Igr^^ , Igj^_ , Igj^^ 
and Ig^^ defined in Eqs. (|38p . p9p . (|43|l . (|44p and (|45p as functions of a and /3. The plots are shown starting from a minimum 
value (left to right) of a„i„ — 0.1, 0.5, 1 up to a maximum value of Qmax = 6. 
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FIG. 5: Integrals of the functions Igx- ' -^axy ' -^sl^ ' -^SLy ^'i*^ 
Ig^^ over the domain ctmin < a < Qmax, —5 < j3 < 5 for 
the different values of a„,i„. We also plot the values of the 
combinations 7z = Igj,^ + Ig^, and ly = Ig^^ + /gj^^^. Notice 
that for all of values of ami„, the largest contribution to 1^ 
comes from Jg^, , that is the wake or transverse mode whereas 
the largest contribution to ly comes from Ig^y, that is the 
sound or longitudinal mode. Moreover, Iz is always larger 
than /„. 



result reads as 



and 




Figure H] shows three dimensional plots of 
^aTz 1 IgTy 1 IgLz 1 IgLy 9.nd Ig^^ as functions of a and 
p. Shown are also the corresponding contour plots. 
To test the sensitivity of the results on the distance 
to the source, the plots are shown starting from a 
minimum value of = 0.1,0.5,1 up to a maximum 
value of a„,ax = 6. Figure [5] shows the integrals of the 
above functions over the domain a„,i„ < a < a^^^^ 
— 5 < /3 < 5 for the different values of a,„i„. Shown also 
are the values of the combinations Iz = IgT- + IgLz 
ly = Ig^ + Igj^ . Notlcc that for all of the shown values 
of the largest contribution to Iz comes from /g^,, 
that is the wake or transverse mode whereas the largest 
contribution to ly comes from Ig^ , that is the sound or 
longitudinal mode. Moreover, Iz is always larger than 
ly. This last result shows that for the case treated here, 
where (^^jv^ < 1, the momentum deposition will be 
preferentially forward peaked. 



iSL)z 



iiir) [dx) V 3r 



dE\ ( 2v 



V47r 



3as 2 
2(7+1) ~ (s + l)2 

2 



dE 
dx 



2v 



Ig.A»,P)- 



(43) 



In a similar fashion we obtain 
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/3 



16 7o '^^Ji^i/^^^^^)^^"*"'^-! 



(s + l)2 2(s + l) s 



(±\ 

Utt/ \ dx ) \?,Ts 



(a,/3). 



(44) 



IV. AZIMUTHAL CORRELATIONS 

We now proceed to use the formalism developed above 
to compute the azimuthal angular correlations for events 
where the leading hadron has a larger or equal mo- 
mentum than the associated ones. Figure [B] shows the 
per-trigger correlations for the cases where the leading 
hadron is produced in a momentum range 3 GcV < pt < 
4 GcV and the associated ones in the momentum ranges 
0.4 GeV < PT < 1 GcV, 1 GeV < pt < 2 GeV, 2 GcV 
< Pt < 3 GeV and 3 GeV < pt < 4 GeV, compared to 
data from PHENIX 0]. To obtain these correlations we 
have generated 2 — >■ 3 hard scattering parton events using 
the MadGraphS event generator [l7| , within the momen- 
tum window 10 GcV < pt < 12 GcV. Out of the three 
partons. we choose the one with the largest momentum 
to become the leading hadron by collinear fragmentation 
using the KKP parton fragmentation functions [ll] . The 
other two partons in the event travel within the medium 
and thus represent the sources of energy-momentum de- 
position in the away-side. We consider that these partons 
travel 6 fm on average [l^ . and a medium induced en- 
ergy loss per unit length dE/dx = 2 GeV/fm. Therefore 
in average, these partons transfer all of their initial mo- 
mentum to the medium. For these away-side partons. 
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FIG. 6: Per-trigger azimuthal angular correlations. Shown 
are four different cases each one having the leading hadron 
produced with a fixed momentum range 3 GeV < pt < 4 
GeV. From top to bottom the associated hadrons are in the 
momentum ranges 0.4 GeV < Pt < 1 GeV, 1 GeV < Pt < 2 
GeV, 2 GeV < pt < 3 GeV and 3 GeV < pt < 4 GeV. The 
results are compared to data from PHENIX Q]. 



Notice that Fig. |6] shows a quite good agreement of 
this simple scenario with the correlation data, particu- 
larly since it reproduces the systematics of double-hump 
decreasing intensity as the momentum of the away-side 
hadrons becomes closer to the momentum of the leading 
hadron. 



V. SUMMARY AND CONCLUSIONS 

In this work, we have studied the way that fast moving 
partons deposit their energy-momentum when traveling 
in a medium. We have shown that for conditions re- 
sembling a medium produced in a heavy ion reaction, 
energy-momentum is preferentially deposited along the 
head-shock and that in order to produce a double-hump 
structure in the away-side of azimuthal correlations, one 
can consider that two partons, instead of one, travel to- 
wards the away-side. We have argued that these away- 
side partons can be originally produced by a 2 — >■ 3 hard 
scattering. Even though 2 — > 3 parton processes arc sup- 
pressed with respect to 2 — > 2 ones by and extra power of 
a^, there is also a kinematic enhancement that partially 
compensates. This is due to the fact that for a given en- 
ergy deposited in the away-side, the partons from 2 — > 3 
processes are produced with lower momentum than the 
away-side parton in 2 — >■ 2 processes, thus the former are 
more copiously produced. 

We have resorted to linear viscous hydrodynamics to 
compute the medium response to the fast moving par- 
tons and the Cooper-Frye formula to compute the par- 
ton distribution originated from this disturbance in a 
static background medium. To mimic a scenario where 



we compute the angular distribution around their direc- 
tion of motion by means of the Cooper-Frye and linear 
viscous hydrodynamics formalisms, described in Sees. HIl 
and mil We use the value Q;„i„ = 0.1. Since we consider 
massless partons, we take their speed v = 0.9995c. The 
sound attenuation length is taken to be = l/(37rTo) 
which results from taking r\ = l/(47r) and consider the 
value To =350 MeV, and the sound velocity Cs = 1/a/3. 
In order to produce final state hadrons out of the away- 
side partons that emerge via the Cooper-Frye formula, we 
consider that the partons in the head shock region carry 
half of the hard scattering parton momentum. We sub- 
sequently let these partons produce final state hadrons 
using once again the KKP in vacuum fragmentation 
functions. This procedure is intended to mimic surface 
emission of the leading parton together with in-medium 
energy-momentum deposited by the away-side ones. The 
hadron momentum distributions obtained from the lead- 
ing and associate partons are shown in Figs. [7] and [51 
respectively. 




4 6 8 
P, [GeV/c] 



FIG. 7: Hadron momentum distribution obtained from a lead- 
ing parton with momentum in the range 10 GeV < pt < 12 
GeV. 
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[GeV/c] 

FIG. 8: Hadron momentum distribution obtained from an 
associated parton with momentum in the range 5 GeV < pr < 
6 GeV. 

the leading hadron conies from surface emission and the 
away-side partons deposit all of their original energy- 



momentum within the medium, we have selected partons 
produced in 2 — >■ 3 processes, evolving the parton with 
the highest energy to become the leading hadron using 
KKP in vacuum fragmentation. Also, the head-shock 
region around the direction of motion of the away-side 
partons, is assumed to carry half of the original hard- 
scattering parton momentum. The partons in the head- 
shock region are then hadronized using also KKP in vac- 
uum fragmentation. The comparison to PHENIX data 
shows that this simple scenario reproduces the system- 
atics of a decreasing away-side correlation when the mo- 
mentum of the associated hadrons becomes closer to the 
momentum of the leading hadron. This scenario seems 
to avoid the shortcomings of the Mach cone as the origin 
of the double-hump structure in the away-side. A more 
detailed analysis is currently under way to explore the 
systematics of away-side correlations for a larger set of 
conditions and will be reported elsewhere. 
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